\documentclass[10pt,a4paper]{article} 

\usepackage{ctex} % 中文支持
\usepackage[top=2.5cm, bottom=2.5cm, left=2.5cm, right=2.5cm]{geometry} % 页边距
\usepackage{amsmath, amssymb} % 数学公式与符号
\usepackage{graphicx}
\usepackage{pythonhighlight}
\usepackage{url} 

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%文档的题目、作者与日期
\author{五六七 }
\title{悬链线的方程 }

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
\begin{document}

\maketitle

\begin{abstract}
建立悬链线符合的微分方程，求解函数表达式。并用观测数据进行验证。
\end{abstract}

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%\setcounter{tocdepth}{2}
%\renewcommand\contentsname{目录}
%
%\renewcommand {\baselinestretch} {1.3}\normalsize 
%\tableofcontents 
%\renewcommand {\baselinestretch} {1.0}\normalsize


%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

\section{问题描述}
一根链子的两段固定在墙上的两个点，其余根据重力自然下垂。求曲线的函数表达式。
用实验得到观测数据，求出函数表达式的参数。

\begin{figure}[ht!]\centering
\includegraphics [height=4cm, width=9cm]{catenary03.png}
\caption{一根真正的悬链线 }
\end{figure}

\begin{figure}[ht!]\centering
\includegraphics [height=4cm, width=9cm]{pantograph-catenary-systemblocks.png}
\caption{悬链线的另一个例子 }
\end{figure}




实验观测：可对实物拍照，然后量出很多点的坐标，然后代入模型函数，求出参数。
或将一张白纸放在墙上，挂一根线，然后描出悬链线的形状，然后同样的数据分析。


%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

\section{建立模型：导出悬链线的微分方程}

\begin{enumerate}\itemsep0.5em

%\item  {\color{red}我们来导出悬链线满足的微分方程、边界条件、并求解。} 

\item  设悬链线的方程为 $y(x)$, 设两端的坐标为 $(x_1,y_1)$ 与 $(x_2,y_2)$.

\item  设悬链线上相邻两点的坐标分别为 $P(x,y(x))$ 和 $Q(x+\Delta x, y(x+\Delta x))$. 

\begin{center}
\includegraphics [height=4cm, width=9cm]{catenary.png}
\end{center}

\item  设$PQ$之间的这段弧长为 $\Delta s$, 其中 $s$ 是从左端点 $P_1$ 到$P$的弧长。

\item  设单位长度悬链线受到的重力为 $\gamma$, 则弧 $PQ$ 受到的重力为 $W=\gamma \Delta s$. 

\item  设在 $P$ $Q$ 两点的张力分别为 $F_1$ 和 $F_2$, 这是两个向量。

\item  设张力在水平方向上的分量分别为 $H_1=H(x)$ 和 $H_2=H(x+\Delta x)$. 

\item  设张力在竖直方向上的分量分别为 $V_1=V(x)$ 和 $V_2=V(x+\Delta x)$. 

\item  根据水平方向受力平衡，可得 $H_1=H_2$, 记为 $H_0$. 

\item  根据竖直方向受力平衡，可得 $V_1+W=V_2$.  由此将导出悬链线方程。

\item  从 $V_2 - V_1 = W$ 可得 $V(x+\Delta x) - V(x) = \gamma \Delta s$. 
两边除以 $\Delta x$ 可得 $$\frac{V(x+\Delta x) - V(x)}{\Delta x} = \gamma \frac{\Delta s}{\Delta x}. $$

\item  根据勾股定理，可得 $\Delta s = \sqrt{(\Delta x)^2 + [y(x+\Delta x)-y(x)]^2}$. 除以 $\Delta x$ 可得
$$\frac{\Delta s}{\Delta x} = \frac{\sqrt{(\Delta x)^2 + [y(x+\Delta x)-y(x)]^2}}{\Delta x}. $$ 

\item  对上述两式求 $\Delta\to 0$ 时的极限，可得 $$V'(x) = \gamma \sqrt{1+[y'(x)]^2}. $$

\item  因为张力的竖直方向与水平方向的比例，正好就是切线的斜率，所以有 $$y'(x) = \frac{V(x)}{H(x)} = \frac{V(x)}{H_0}.$$

\item  从上述两式可得悬链线满足的微分方程 $$ y'' = \frac{V'(x)}{H_0} = \frac{\gamma}{H_0}\sqrt{1+(y')^2}. $$

\item  记 $\frac{\gamma}{H_0}=a$. 可将方程简写成 $$y'' = a\sqrt{1+(y')^2}. $$

\end{enumerate}

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

\section{求解模型：求解悬链线方程}

\begin{enumerate}\itemsep0.5em

\item  记 $z=y'$, 所求方程化为 $z' = a\sqrt{1+z^2}$. 

\item  分离变量可得 $\frac{dz}{\sqrt{1+z^2}} = adx$. 两边积分可得 %$$\int \frac{dz}{\sqrt{1+z^2}} = 
$$\text{arcsinh} z= ax +C_1. $$

\item  由此可得 $z$ 是 $x$ 的双曲正弦函数 $z=\text{sinh}(ax+C_1)$. 

\item  积分可得 $y$ 是 $x$ 的双曲余弦函数 $y=\frac{1}{a}\cosh (ax+C_1) + C_2$.

\item  其中双曲正弦函数与双曲余弦函数分别为 
$$\text{sinh}(x) = \frac{e^x - e^{-x}}{2},\,\,\, \text{cosh}(x) = \frac{e^x + e^{-x}}{2}. $$ 

\end{enumerate}

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

\section{求解模型：边值条件}

\begin{enumerate}\itemsep0.5em

\item  因为悬链线的两端是固定的，所以有边值条件 
\begin{eqnarray*}
\left\{\begin{array}{rcl}
y(x_1) &=& y_1, \\ 
y(x_2) &=& y_2. 
 \end{array}\right. 
\end{eqnarray*}

\item  将边值条件代入解函数，可得
\begin{eqnarray*}
\left\{\begin{array}{rcl}
\cosh (ax_1 + C_1) + C_2 &=& ay_1, \\
\cosh (ax_2 + C_1) + C_2 &=& ay_2. 
 \end{array}\right. 
\end{eqnarray*}

\end{enumerate}

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

\section{求解模型：从悬链线方程确定悬链线方程的参数}

\begin{enumerate}\itemsep0.5em

\item  设悬链线的长度为 $L$, 则有 $L>\sqrt{(x_2-x_1)^2+(y_2-y_1)^2}$.

\item  根据弧长公式，可得 $$L=\int_{x_1}^{x_2} \sqrt{1+[y'(x)]^2}dx. $$

\item  根据悬链线方程 $y'' = a\sqrt{1+(y')^2}$, 可得上述积分为 $$L = \frac{1}{a}\int_{x_1}^{x_2} y''(x)dx = \frac{1}{a}[y'(x_2)-y'(x_1)].$$

\item  代入 $y'=z=\text{sinh}(ax+C_1)$ 可得 $$L=\frac{1}{a}[\text{sinh}(ax_2+C_1) - \text{sinh}(ax_1+C_1)]. $$

\item  根据双曲正弦的和差化积公式，可得
$$ aL = 2\text{sinh}\left[ \frac{a(x_2-x_1)}{2} \right] \text{cosh} \left[ \frac{a(x_2+x_1)+2C_1)}{2} \right]. $$

\item  将两个边值条件相减，并由双曲余弦的和差化积公式，可得
\begin{eqnarray*}
a(y_2-y_1) &=& \text{cosh}(ax_2+C_1) - \text{cosh}(ax_1+C_1) \\ 
&=& 2\text{sinh}\left[ \frac{a(x_2-x_1)}{2} \right] \text{sinh} \left[ \frac{a(x_2+x_1)+2C_1)}{2} \right].
\end{eqnarray*}

\item  根据 $\text{cosh}^2(x) - \text{sinh}^2(x) = 1$ 可得 
$$\frac{a}{2}\sqrt{L^2 - (y_2-y_1)^2} = \text{sinh}^2 \left[ \frac{a(x_2-x_1)}{2} \right]. $$


\end{enumerate}

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

\section{求解模型：从实际观测数据确定参数}

\begin{enumerate}\itemsep0.5em 

\item   从微分方程模型可得悬链线的函数关系式是双曲余弦函数 
$$y=\frac{1}{a}\cosh (ax+C_1) + C_2.$$

\item  观测得到一些坐标数据 $(x_1,y_1), (x_2,y_2),\cdots, (x_n,y_n)$. 

\item  使用最小二乘法求出参数 $a,C_1,C_2$ 的值。

\end{enumerate}

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

\section{编程计算}





%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

\section{回答问题}


%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%\section{习题}
%导出悬链线满足的微分方程，并写出边界条件。



%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%\section{参考文献 }
\begin{thebibliography}{99}
\bibitem{dingtongren} 丁同仁、李承治，\emph{常微分方程教程}，高等教育出版社，2022年3月第三版。
\bibitem{sishoukui-2} 司守奎,孙玺菁. \emph{Python数学建模算法与应用}, 国防工业出版社. 2022年1月第1版. 


\end{thebibliography}

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

\end{document}

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%





